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*  * 


1  Purpose  of  the  Work 

The  overall  goal  of  this  project  is  to  develop  an  accurate  design  tool  for  predicting  and  con¬ 
trolling  unsteady  oscillations  in  high-performance,  gas-turbine  combustors.  The  sensitivity 
equation  method  (SEM)  has  been  developed  in  a  stand-alone  commercial  package  (called 
SENSE)  to  the  point  of  providing  turbulent  flow  sensitivities  for  chemically  reacting  flows 
(and  all  subsets).  This  STTR  project  adds  the  element  of  unsteadiness  which  no  sensitivity 
method  exploits  better  than  SEM.  Our  focus  during  Phase  I  is  to  develop  such  a  tool  using 
the  GASV  and  SENSE  CFD  software  to  study  the  aero-acoustic  instabilities  observed  in 
the  National  Energy  Technology  Laboratory  (NETL)  combustor.  Ultimately  this  software 
will  be  useful  for  determining  optimal  placement  of  controllers  and  transducers. 

1.1  Plans  for  Phase  II 

After  the  Phase  I,  we  propose  to  begin  a  three-part  effort  in  Phase  II  -  each  directed  by  a 
member  of  the  STTR  team.  AeroSoft’s  portion  will  focus  on  developing  a  driver  for  coupling 
both  the  flow  and  sensitivity  solvers.  Much  of  the  computational  work  performed  in  a  flow 
solver  can  be  re-used  to  compute  sensitivities.  These  include  the  flow  solution  (of  course) 
as  well  as  grid  area,  volume  and  metrics,  flux  Jacobians,  time-integration  methods.  For 
unsteady  applications,  the  driver  would  be  used  to  call  GASViox  a  flow  update  and  then 
pass  required  information  for  a  sensitivity  update.  In  the  steady  case,  the  driver  would 
call  QASV  once  the  flow  solution  has  converged  and  then  pass  the  relevant  information 
to  SENSE.  The  majority  of  this  exchange  is  controlled  using  a  Message  Passing  Interface 
(MPI).  This  solver  coupling  increases  the  marketability  of  the  sensitivity  software  because 
the  connection  between  flow  and  sensitivity  solver  is  completed  automatically.  Currently, 
the  user  must  be  able  to  export  a  specific  solution  format  and  then  set  up  input  files  for 
SENSE.  The  coupling  software  therefore  simplifies  the  user’s  workload.  If  a  generic  flow 
solver  is  selected,  then  software  currently  existing  in  SENSE  would  be  used  to  compute  the 
re-usable  information. 

The  ICAM  effort  would  focus  on  using  the  Ansatz  from  Floquet  theory  to  construct  a 
solution  to  the  linear  sensitivity  equations  for  a  periodic  flow.  We  would  exploit  this  structure 
for  efficient  computation  of  this  class  of  time-varying  flow  sensitivity  problems.  Also,  ICAM 
will  work  on  sensitivity  problems  involving  turbulent  flows  where  the  design  variable  controls 
the  geometry.  This  is  the  biggest  stumbling  block  in  SENSE  currently  and  would  greatly 
increase  the  marketability  of  the  product. 

The  Mechanical  Engineering  Department  would  be  responsible  for  designing  and  com¬ 
pleting  an  unsteady,  combustion  experiment  of  broad  interest  which  could  be  simulated  with 
the  current  capabilities  of  QASV .  Rather  than  request  funding  to  increase  the  fidelity  of 
QASV  to  run  a  specific  combustion  problem,  we  feel  a  more  productive  approach  is  to  design 
a  combustion  problem  where  the  focus  can  be  placed  on  the  sensitivity  problem.  The  exper¬ 
imental  geometry  and  results  would  serve  as  a  test-bed  for  the  driver/boundary-condition 
software  developments  given  above.  Also,  a  team  member  will  be  responsible  for  running 
an  unsteady  QASV/SENSS  simulation  in  the  current  uncoupled  fashion  to  compare  with 
the  improvements  provided  by  the  coupling  driver  software.  This  comparison  would  occur 
at  the  end  of  the  Phase  II. 
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1.2  Elements  of  the  Phase  I 

During  the  12  months  of  the  project,  the  following  technical  objectives  have  been  accom¬ 
plished: 

1.  We  implemented  a  simple,  yet  adequate  chemistry  mechanism  into  QASV  and  SENSE, 
and  applied  to  the  reacting  flow  within  the  NETL  rig.  With  a  quasi-steady  flow 
solution,  we  computed  the  sensitivity  to  inlet  velocity  and  swirler-to-injector  distance. 

2.  Using  the  quasi-steady  solution  as  a  starting  point,  we  simulated  the  unsteady  com¬ 
bustion  within  the  NETL  rig.  We  determined  the  requirements  to  accurately  simulate 
the  unsteady  oscillations  induced  by  thermo-acoustic  interaction. 

3.  Using  time-frozen  snapshots  from  the  unsteady  flow  simulation,  we  computed  sen¬ 
sitivities  to  inlet  velocity  and  swirler  location.  Note  that  the  Sensitivity-Equation 
Method  (SEM)  yields  a  partial  differential  equation  for  the  unsteady  sensitivity.  Our 
goal  in  this  STTR  project  is  to  utilize  this  capability,  which  other  sensitivity  methods 
lack. 

4.  We  solved  the  one-dimensional,  unsteady  flow  and  sensitivity  equations  simultaneously. 
The  initial  value  problem  corresponds  to  a  10  :  1  (by  pressure)  Riemann  problem.  Also, 
a  quasi- ID  approximation  is  implemented  to  determine  the  sensitivity  to  inlet  density 
and  geometric  changes. 

1.3  Conclusions  from  the  Phase  I 

Several  levels  of  sophistication  in  chemistry  modeling  have  been  investigated  including  one/two- 
step  models,  reduced  global  models  and  fundamental  mechanisms.  We  conclude  that  a  model 
which  uses  a  fundamental  reaction  mechanism  provides  the  best  performance  both  numeri¬ 
cally  and  physically.  We  also  have  determined  the  best  way  to  solve  an  internal  combustion 
problem  where  no  geometric  flame  holder  exists,  such  as  the  NETL  combustor.  We  have 
computed  flow  sensitivities  to  the  inlet  velocity  and  swirler-ring  location.  These  sensitivities 
are  most  significant  around  the  flame  front.  The  feasibility  of  a  time-accurate  sensitivity 
solution  is  pro/en  using  a  one-dimensional  Riemann  problem. 
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2  Brief  Description  of  the  Work  Carried  Out 

Our  Phase  I  goal  is  to  simulate  a  forced  instability  in  a  simplified  version  of  the  NETL 
combustor.  A  baseline  pre-mixed  profile  is  imposed  as  a  planar  jet  into  the  combustion  region. 
This  profile  is  determined  by  solving  the  steady-state,  axi-symmetric  flow  equations  for  the 
fuel  nozzle  alone.  To  determine  the  sensitivity  profile,  we  assume  the  mass-fraction  profiles 
can  be  approximated  using  a  cubic  Lagrange  polynomial  in  the  swirler  location.  Turbulent 
mixing  of  air  and  methane  in  the  nozzle  occurs  via  a  two-equation  model.  The  baseline 
profile  is  made  unsteady  via  a  time-dependent  sinusoidal  function  where  the  amplitude  and 
frequency  are  deduced  from  experimental  data. 

The  flow  in  the  combustor  uses  a  finite-rate  approximation  with  twelve  reacting  species. 
The  steady- state  solution  serves  as  the  starting  point  (i.e.,  t  =  0)  for  the  time-dependent 
calculation.  At  several  points  in  time,  the  flow  variables  are  output  and  a  sensitivity  calcula¬ 
tion  is  performed.  While  this  is  not  a  true  unsteady  sensitivity,  we  have  established  the  best 
chemistry  model,  included  the  most  important  physics  in  the  study  and  gained  sensitivity  in¬ 
sight.  A  true  unsteady  sensitivity  solution  has  been  completed  by  solving  a  one-dimensional 
Riemann  problem  and  a  quasi-ID  problem. 

3  Results 

3.1  The  Design  Problem 

Four  different  swirler  configurations  are  available  in  the  fuel-nozzle  section  of  the  NETL 
gas-turbine  combustor.  In  particular,  the  swirl  vanes  can  be  placed  farther  upstream  of  the 
fuel-spoke  injector  in  increments  as  shown  in  Fig.  1.  The  first  case  corresponds  to  locating 
the  swirl  vanes  3.25”  upstream  of  the  combustor.  Each  successive  case  corresponds  to  the 
vanes  being  located  one  inch  farther  upstream  (i.e.,  to  the  left).  One  of  our  design  variables 
is  this  swirler  location. 

Under  the  calendar  and  financial  constraints  of  a  Phase  I  effort,  we  seek  to  simulate  a 
chemically  reacting  flow  which  captures  the  important  physics  of  the  gas-turbine  combustor 
while  excluding  more  complicated  features  such  as  a  self-excited  flow  instability  and  three- 
dimensional  effects.  Our  goal  is  not  lo  model  the  entire  NETL  rig,  but  to  identify  the 
relevant  aspects  of  the  flow  and  then,  investigate  the  sensitivity  to  the  swirler’s  location.  We 
believe  that  these  important  features  include  (1)  swirling  fuel/air  mixing  in  the  fuel  nozzle, 
(2)  turbulent  combustion  of  a  pre-mixed  flame  in  the  combustor,  and  (3)  induced  unsteady 
mass  flow  which  leads  to  an  unsteady  fuel-to-air  equivalence  ratio. 

For  computational  simplicity,  we  approximate  the  three-dimensional  full  NETL  config¬ 
uration  with  a  two-dimensional  version  which  begins  at  the  nozzle  exit  and  extends  to  the 
end  of  the  refractory  plug.  After  trying  several  simplified  methane-air  chemistry  models 
including  one/two-step  mechanisms  [1]  and  reduced  global  models  [2,  3],  we  resort  to  a 
more  fundamental  reaction  mechanism  [4]  which  contains  12  species  and  32  elementary  re¬ 
actions  (including  collision- partner  reactions  separately).  After  testing,  the  one/ two-step 
mechanisms  performed  with  extreme  stiffness.  Increasingly  small  time  steps  failed  to  yield 
a  reasonable  solution.  The  complexity  of  the  Paczko/Peters  reduced  models  justified  using 
numerical  source  Jacobians  rather  than  tediously  deriving  analytical  ones.  This  added  to  the 
computational  cost  per  iteration  and  required  major  modifications  to  the  sensitivity  software. 
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Figure  1:  Fuel  nozzle  showing  the  range  of  positions  for  the  swirling  vanes. 


In  addition,  the  reduced  models  required  a  small  trace  amount' of  product  concentration  in 
order  for  the  reactions  to  occur.  Imposing  a  minimum  floor  upon  the  product  species  further 
hampered  the  models’  computational  effectiveness. 

In  contrast,  the  format  of  the  twelve-specie  model  fit  well  into  the  structure  previously 
developed  in  AeroSoft’s  flow  and  sensitivity  solvers.  Therefore,  we  eventually  abandoned 
the  reduced  global  models  in  favor  of  a  model  based  on  a  fundamental  reaction  mechanism. 
The  drawback  of  this  approach  is  that  the  larger  number  of  species  significantly  added  to 
the  computational  cost. 

3.2  Swirler/Spoke  Nozzle  Flow 

To  determine  the  flow  profile  at  the  inlet  to  the  combustion  region,  we  simulate  the  three- 
dimensional  mixing  of  air  and  fuel  that  takes  place  in  the  fuel  nozzle.  This  is  a  separate 
up-front  calculation  where  the  flow  is  assumed  frozen  and  turbulent  (a  two-equation  K—u> 
model  is  used).  The  mass-fraction  profiles  for  N%,  O2  and  CH4  are  then  applied  as  an  in-flow 
profile  for  the  two-dimensional,  chemically  reacting  simulation  in  the  combustor.  By  running 
all  four  swirler-ring  cases,  we  can  formulate  the  sensitivity  of  the  mass-fraction  profiles  to 
the  swirler  location  through  a  Lagrange  polynomial.  Along  with  the  flow  variables,  the 
sensitivity  equations  can  be  solved  given  this  mass-fraction  sensitivity  profile. 

Prateep  Chatterjee,  a  doctoral  student  in  the  Mechanical  Engineering  Department  at 
Virginia  Tech  under  the  direction  of  Dr.  Uri  Vandsburger,  has  computed  the  mixing  which 
occurs  in  the  fuel  nozzle.  Air  is  introduced  at  a  swirl  angle  of  45°  and  mixes  with  fuel  at 
the  spoke  ring.  Internal  mass  sources  tuned  for  an  equivalence  ratio  of  <j>  =  0.74  are  used  to 
introduce  the  fuel  into  the  stream.  The  mass-fraction  profiles  and  axial  velocity  at  the  exit 
of  the  fuel  nozzle  are  shown  in  Fig.  2.  As  the  swirling  rings  are  located  farther  upstream  of 
the  combustion  region,  the  swirl  ratio  decreases  at  the  spoke-ring  location.  As  a  result,'  the 
mixing  in  Case  4  is  less  than  in  Case  1. 
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Normalized  distance,  y/h 


(a)  Methane  mass-fraction  profiles.  (b)  Oxygen  mass-fraction  profiles. 


(c)  Axial  velocity  profiles. 


Figure  2:  Profiles  from  the  axi-symmetric  calculations  where  the  swirler  location  progresses 
upstream  in  each  of  the  four  cases. 


Given  the  four  mass-fraction  profiles  for  specie  i  at  the  nozzle  exit,  the  dependency  of 
the  local  mass  fraction  upon  the  swirler-spoke  distance  (r/s)  can  be  expressed  as 


Yi(y;ris)  =  'ZYi 


n  vs -vs 

k= i 

U) 


(k) 


n 

k=l 


(i) 


The  above  mass-fraction  profile  is  a  cubic  function  of  the  swirler  distance.  Carpet  plots  of 
Ycha  and  Yo2  and  their  sensitivities  to  r/s  are  given  in  Fig.  3.  The  sensitivity  of  the  profile 
can  be  determined  by  differentiation  of  the  Lagrange  polynomial  with  respect  to  r/s.  The 
trend  discussed  above  relating  decreasing  mixing  with  increasing  distance  is  evident  by  the 
negative  sensitivity  near  y  =  0  and  y  =  1  for  r)s  =  0-75  (see  Fig.  3(b)).  That  is,  if  the 
swirler-spoke  distance  increases,  then  the  mass  fraction  of  CH4  will  decrease  at  the  fringes 
of  the  mixing  layer  meaning  the  fuel  has  mixed  less. 

We  note  that  early  in  the  project  we  proposed  to  introduce  methane  via  a  normal- 
injection  hole.  We  could  then  determine  the  sensitivity  of  the  flow  to  the  hole  location  by 
differentiating  a  continuous  flow  profile.  This  idea,  while  abandoned  in  the  Phase  I  because 
it  leads  to  a  diffusion  flame  in  the  combustor  rather  than  a  pre-mixed  flame,  has  produced 
an  AIAA  paper  which  will  be  presented  in  January,  2002  [5]. 


3.3  Combustion-Region  Flow 
3.3.1  Solution  Description 

Using  AeroSoft’s  QASV  flow  solver,  the  NETL  flow  is  simulated  with  third-order  spatial  ac¬ 
curacy  and  second-order  temporal  accuracy.  Temporal  accuracy  is  obtained  using  a  dual-time 
stepping  procedure.  To  take  reasonably  large  physical  time  steps,  the  stiff  specie  production 
terms  (i.e.,  large  in  magnitude)  must  be  treated  implicitly.  Using  explicit  schemes  to  solve 
numerical  problems  with  large  source  terms  requires  a  very  small  time  step  to  maintain  sta¬ 
bility.  The  implicit  procedure  in  the  dual-time  stepping  method  allows  us  to  obtain  answers 
in  much  fewer  time  steps.  The  caveat  is  that  an  inner-iteration  problem  must  be  converged 
for  each  physical  time  step. 

The  chemistry  model  utilizes  the  methane-air  mechanism  described  in  Bowman  and 
Seery’s  article  [4]  on  shock-induced  combustion.  The  data  for  the  model  is  given  in  Ta¬ 
ble  1 .  The  thermodynamic  properties  for  each  of  the  twelve  species  are  computed  using  the 
LeRC  curve  fits  compiled  by  McBride  [6,  7].  By  casting  the  chemical  equilibrium  problem  as 
a  minimization  of  Gibbs  free  energy,  we  can  use  McBride’s  data  to  determine  the  equilibrium 
constant  and  thus,  the  backward  rate. 

A  number  of  ideas  were  utilized  to  advance  the  simulation  to  the  steady  state.  First,  we 
mention  two  ideas  that  were  not  effective. 

1.  The  inlet  temperature  of  the  methane/air  mixture  is  too  low  (T)  =  583  K)  to  ignite. 
Furthermore,  sources  for  combustion  do  not  exist  anywhere  in  the  flow  field  (e.g. ,  a 
hot  wall  or  shock-induced  combustion).  A  heat  source  is  necessary  to  ignite  the  flow, 
yet  none  exists  physically.  The  result  is  a  cold-flow  calculation. 
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(c)  Oxygen  mass  fraction,  Yo2  {y,  Vs)-  (d)  Oxygen  sensitivity,  dYo2 /dr), . 

Figure  3:  The  mass  fraction  profiles  and  sensitivities  in  (t/,r/s)  space. 
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Reaction 

kf  =  AT*1  exp(- 

EJT ) 

3 

A 

V 

Ea 

1 

CHi  +  M  #  CH3  +  H  +  M 

2.0  x  1017 

0.0 

44, 500 

2 

CHa  +  OHr±  CHz  +  H20 

2.8  x  1013 

0.0 

2,500 

3 

CH±  +  CH%  +  OH 

2.0  x  1013 

0.0 

4,640 

4 

CHa  +  H  ^  CHz  +  H2 

6.9  x  1013 

0.0 

5, 950 

5 

CHZ  +  02^  HCO  +  H20 

2.0  x  1010 

0.0 

0 

6 

CHz  +  0^  HCO  +  H2 

1.0  x  1014 

0.0 

0 

7 

HCO  +  OH  #  CO  +  H20 

1.0  x  1014 

0.0 

0 

8 

HCO  +  M  ^  H->- C  O +  M 

2.0  x  1012 

0.5 

14, 400 

9 

CO  +  OH  ^  C02  +  H 

5.6  x  1011 

0.0 

545 

10 

H  +  02^0  +  0H 

2.2  x  1014 

0.0 

8, 340 

11 

o  +  h2^  h  +  oh 

1.7  x  1013 

0.0 

4, 750 

12 

O  +  H20  ^  2 OH 

5.8  x  1013 

0.0 

9,070 

13 

H  +  H20  ^  H2  +  OH 

8.4  x  1013 

0.0 

10, 100 

14 

H  +  OH  +  M^  H20  +  M 

1.0  x  1019 

-1.0 

0 

Table  1:  Methane-Oxygen  reaction  mechanism  utilized  in  the  NETL  combustion  region 
(from  Bowman  and  Seery  [4]).  Collision  efficiencies  cited  from  Warnatz  [8]  are  M  = 
H2  [1.0],  02  [0A},N2  [0A},H2O  [6.5],  CO  [0.75],  C02  [1.5],  CH4  [6.5],  Note:  N2  is  an  in- 
ert  specie. 

2.  Raising  the  static  temperature  artificially  by  manipulating  the  specie  densities  at  cer¬ 
tain  strategic  spots  in  the  flow  results  in  local  “explosions”.  The  local  velocity  pattern 
is  corrupted  significantly  and  the  local  Mach  number  increases  unrealistically.  After 
disabling  the  temperature  source,  the  flow  field  does  not  regain  a  realizable  pattern 
and  the  approach  is  abandoned. 

To  obtain  the  correct  flow  pattern^e  must  work  backwards  beginning  with  the  primary 
products  rather  than  the  reactants.  The  flow  is  assumed  to  consist  only  of  N2,  C02 ,  and 
H20  at  a  high  temperature,  T  =  1500  K.  Once  this  frozen  flow  calculation  has  converged, 
we  convert  the  solution  to  the  12-specie  finite-rate  chemistry  model.  We  then  impose  the 
inlet  conditions  as  a  boundary  condition  using  the  hot-flow  solution  as  an  initial  restart. 

The  inlet  conditions  for  the  nozzle  correspond  to  data  provided  by  Douglas  Straub  of 
the  National  Energy  Technology  Laboratory.  The  equivalence  ratio,  nozzle  velocity,  and 
inlet  temperature  are  <j)  =  0.74,  U{  =  55.6  m/s  and  T{  =  583  Ah  The  total  pressure  in  the 
plenum  is  p0  =  7.5  atm,  and  this  results  in  an  inlet  Mach  number  of  M,  =  0.1172.  Once 
the  flow  with  these  uniform  inlet  conditions  becomes  pseudo-steady,  we  apply  the  profile 
from  the  axi-symmetric  mixing  calculations  discussed  in  Sec.  3.2.  Flow  visualization  of  the 
temperature,  Mach  number,  and  reactants  are  shown  in  Fig.  4. 

The  flow  consists  of  a  planar  pre-mixed  jet  which  generates  two  counter-rotating  vortices 
-  one  at  the  base  region  of  the  combustor  and  the  other  upstream  of  the  constrictor  plug. 
The  vortices  squeeze  the  flame  which  is  completely  burned  approximately  2.6”  into  the 
combustor. 
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(c)  Methane  mass  fraction,  Ych4- 


(d)  Oxygen  mass  fraction, 


Figure  4:  Flow  visualization  of  the  steady  flow  solution. 
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3.3.2  Integral  Sensitivity  Equations 

The  integral  form  of  the  time-dependent  sensitivity  equations  can  be  written  as  follows. 

JJJ  S  dV  +  J  F'(S;  Q)  •  ndA  —  J^  F(,(S;  Q)  •  ncL4  +  JJJ  W'  dV  (2) 

The  flow  sensitivity  is  denoted  as  S  and  is  defined  as  dQ/dr).  Q  represents  the  vector 
of  conserved  variables  and  77  is  a  generic  design  variable.  The  sensitivity  of  the  inviscid 
and  viscous  flux  vectors  are  denoted  as  F7  and  F^,  respectively.  The  change  of  chenv.cal 
composition  by  combustion  is  controlled  by  the  source  term,  W.  Before  the  sensitivities 
can  be  determined,  the  flow  variables  must  be  known  in  the  form  of  Q.  In  our  “snap-shot” 
sensitivities,  we  solve  the  steady-state  version  of  Eqn.  (2). 

3.4  Unsteady  Flow  History 

In  the  NETL  experiments,  a  lean  mixture  of  fuel/air  burns  on  the  edge  of  being  extinguished 
in  the  combustor.  As  a  result,  pressure  waves  oscillate  from  the  combustion  region  and  affect 
the  mass  flow  of  air  (rhair)  entering  the  nozzle  from  the  plenum.  (Note:  the  mass  flow  of 
fuel  is  choked  and,  thus,  fixed  in  time.)  The  changing  mass  flow  generates  a  time-dependent 
equivalence  ratio  at  the  flame  front.  As  this  mixture  burns,  the  pressure  waves  increase 
in  amplitude  until  the  oscillations  resonate  at  the  vibrational  modes  of  the  system.  If  this 
process  continues,  the  combustor  can  be  damaged.  The  key  component  that  we  focus  upon 
is  that  the  aero-acoustic  instability  causes  the  mass  flow  of  air  to  be  unsteady. 

In  our  investigation,  we  separate  the  effects  of  the  aero-acoustic  instability  by  analyzing 
a  forced  unsteadiness  at  the  inlet  to  the  combustor.  In  particular,  by  making  the  mass 
fraction  profiles  a  function  of  time  we  can  vary  the  equivalence  ratio  without  affecting  mass 
flow.  Likewise,  we  can  change  the  flow  velocity  without  changing  the  equivalence  ratio. 
Finally,  we  can  change  the  density  profiles  and  change  both  simultaneously.  By  denoting  the 
baseline  properties  with  a  0  superscript,  we  can  summarize  the  cases  as  follows. 

•  Variable  mass  flow,  Fixed  equivalence  ratio. 

PN2  =  p°V)v2(y;^) 

po2  =  P°YS2(y;vs) 

Pch2  =  p°YcH2(y^s) 

u  =  u°(y,  T]s)  [1  -I-  esin(wf)] 

•  Fixed  mass  flow,  Variable  equivalence  ratio. 

pn2  =  p°>n2(i/;^)[1  +  esinM)] 

po2  =  P°Yo2{y,Vs)[i  +  esm(ujt)} 

Pch2  =  P°  -  PN2  -  po2 
u  =  u°(y\ T]s) 
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•  Variable  mass  flow,  Variable  equivalence  ratio. 

Pn2  =  p°Y$2(y;ris)[l  +  esin(ujt)} 

Po2  =  P°>o2(y;7?s)[l  +  esin(a;f)] 

PC  Hi  =  P°YcH2(y'’Tls) 
u  =  tt°(y;  T)s) 

In  our  case,  the  baseline  mixture  density  is  constant  for  all  {y,Tjs)  combinations  and  equals 
pO  _  4.761  kg/m?.  For  an  equivalence  ratio  of  4>  —  0.74,  the  mass  fractions  are  Yn2  = 
0.7353,  Yo2  =  0.2233  and  YCha  =  0.0414.  The  sinusoidal  oscillation  is  specified  according 
to  the  experimentally  observed  fluctuations  in  pressure.  A  6.93%  change  in  the  amplitude 
corresponds  to  e  =  0.0693  and  the  measured  frequency  is  u  =  237  Hz.  Four  snap-shots  of  the 
temperature  are  shown  in  Fig.  5.  The  methane  mass-fraction  contours  are  given  in  Fig.  6. 
The  flame  bends  toward  the  centerline  over  time,  but  at  a  much  lower  frequency  than  the 
forced  inflow  oscillation. 

3.4.1  Snapshot  Sensitivities 

For  the  three  forced  instabilities  discussed  earlier,  the  sensitivity  profile  is  determined  by 
differentiation  with  respect  to  r?s.  The  resulting  profiles  are  given  below. 

•  Variable  mass  flow,  Fixed  equivalence  ratio. 

p'n,  =  p^Y^nv-.i.) 

Po ,  =  P° 
p'c„,  =  P°  [Yga 
u'  =  0 

•  Fixed  mass  flow,  Variable  equivalence  ratio. 

p'n2  =  p°[^v2]'(y;%)[1  +  €sinM] 
p'o2  =  p°[>o2]'(y;^)[1  +  esinM)] 

Pc  Hi  —  ~p'n2  -  p'o2 

v!  =  0 

•  Variable  mass  flow,  Variable  equivalence  ratio. 

Pn2  =  P°[^iv2],(y;^)[l  +  €sin(cut)] 
p'o2  =  p°[yo2],(y;7?5)[1  +  €sinM)l 

p'c Hi  —  p°  [YcHil'iy'i  Vs) 
u  =  0 

During  the  Phase  I,  we  investigate  only  the  first  case  where  the  mass  flow  is  variable  but  the 
equivalence  ratio  remains  constant. 
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(a)  Time,  t  =  0.0000  sec.  (b)  Time,  t  -  0.0494  sec. 


Time  (sec) 

(e)  Boundary  condition  history. 
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Figure  5:  Unsteady  temperature  contours,  T. 


1.  Sensitivity  to  inlet  velocity,  77  =  U{. 

Fig.  7  shows  the  snap-shot  sensitivities  to  the  inlet  velocity.  One  would  reason  for 
a  fixed  flame  speed,  increasing  the  inlet  velocity  would  move  the  flame  front  farther 
into  the- combustion  region.  This  would  be  indicated  by  a  positive  sensitivity  in  the 
methane  mass  fraction  aft  of  the  flame  front.  Fig.  7(a)  indicates  a  trend  which  moves 
the  flame. 

2.  Sensitivity  to  swirler  location,  77  =  rjs. 

Fig.  8  depicts  the  mass-fraction  sensitivity  to  increasing  uie  swirler  distance  (i.e.,  mov¬ 
ing  the  swirler  farther  upstream  of  the  combustion  region,)-.  For  all  time  levels,  we  see 
the  effects  of  decreased  mixing  in  the  sensitivity  plots.  We  see  that  methane  mass- 
fraction  at  the  core  of  the  mixture  increases  with  increasing  distance.  The  effect  is 
more  pronounced  at  the  last  time  level  studied  than  at  the  first. 

Fig.  9  shows  pressure-sensitivity  contours  with  increasing  swirler  distance,  from  these 
plots  we  notice  that  “pockets”  of  pressure  fluctuation  exist  near  the  flame  front  as  the 
swirler  is  moved.  The  pressure  away  from  the  flame  seems  to  be  relatively  unaffected 
by  a  new  swirler  location. 
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Forced  Fluctuation 


(c)  Time,  t  =  0.0888  sec. 


(d)  Time,  t  =  0.1480  sec. 


u(tyu0)  -  l +es\r\(m) 


Time  (sec) 


(e)  Boundary  condition  history. 

Figure  7:  Contours  of  unsteady  methane  mass- fraction  sensitivity  to  inlet  velocity,  V 


(a)  Time,  t  —  0.0000  sec. 


(b)  Time,  t  =  0.0494  sec. 


(c)  Time,  t  =  0.0888  sec. 


(d)  Time,  t  =  0.1480  sec. 


u(tyu0)  =  ;+£Sin(cor) 


(e)  Boundary  condition  history. 


Figure  8:  Contours  of  unsteady  methane  mass-fraction  sensitivity  to  swirler  location, 

Yn u .  i  V  =  Vs- 
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(a)  Time,  t  =  0.0000  sec. 


(b)  Time,  t  =  0.0494  sec . 


(c)  Time,  t  =  0.0888  sec. 


(d)  Time,  t  =  0.1480  sec. 


u(tyu0i  -  y+esinftorj 


(e)  Boundary  condition  history. 

Figure  9:  Contours  of  unsteady  pressure  sensitivity  to  swirler  location,  p\  17 


3.5  1-D,  Unsteady  Sensitivity  Solutions 

The  above  snap-shot  sensitivities  do  not  incorporate  sensitivity  history  like  a  true  time- 
dependent  sensitivity  solution.  To  compute  real  unsteady  sensitivities,  we  must  solve  the 
sensitivity  equations  immediately  after  taking  a  physical  time  step  in  the  flow.  In  other 
words,  we  must  solve  the  flow/sensitivity  equations  as  a  coupled  set.  We  test  the  coupled 
solver  by  solving  several  test  problems.  The  first  is  an  example  which  compares  well  with 
previously  published  results  by  Appel  and  Gunzburger  [9] .  The  remaining  examples  demon¬ 
strate  sensitivity  with  respect  to  (a)  the  initial  density  (with  a  uniform  heat  source),  (b)  a 
Bezier  parameter  which  describes  the  shape  of  the  nozzle  cross-sectional  area. 


3.5.1  10:1  Riemann  Problem 


To  verify  our  solution  procedure,  we  solve  a  Riemann  problem  proposed  by  Sod  [10].  The 
unsteady  sensitivity  to  the  high-pressure  is  computed  simultaneously.  For  this  example,  we 
consider  a  constant  area  duct  of  unit  length  x  €  (0, 1)  with  initial  conditions  stated  in  terms 
of  density,  velocity  and  pressure  as 


P 

u 

P 


M 

>  =< 

f  0.125  ) 
° 

l  p  J 

'  R 

l  0.1  J 

(3) 


The  initial  flow  is  separated  by  a  “diaphragm”  located  at  x  —  0.5,  and  the  initial  value 
problem  solved  with  100  cells  of  uniform  width.  All  initial  sensitivities  are  uniformly  zero, 
except  pL  =  1.0.  Results  at  t  =  0.142  with  a  time  step  of  At  =  10"3  are  given  in  Fig.  10 
and  Fig.  11.  Note  that  the  sensitivities  compare  well  with  central  difference  solutions  with  a 
step  size  of  A p  =  0.001  (see  Figure  12).  The  noticeable  discrepancies  occur  at  the  locations 
of  singularities  in  the  true  solution. 


3.5.2  1-D  Nozzle  Sensitivities 

In  this  problem,  we  consider  a  linear^  expanding  nozzle  with  constant  initial  data  as.  given 
below. 

p(0,  x)  =  1.2256,  pu{ 0,  x)  =  833.65,  pe0(0,  x)  =  5.3667  x  105 

We  consider  the  sensitivity  of  the  evolving  flow  to  the  inlet  density  with  a  uniform  heat 
source.  The  idea  behind  the  heat  source  is  to  generalize  the  specie  production  in  the  1-D 
case.  The  results  at  t  =  0.142  sec  are  given  in  Fig.  13.  If  we  specify  the  shape  of  the  nozzle 
using  a  Bezier  curve,  we  can  determine  the  sensitivity  to  a  control  point.  The  results  for  this 
problem  are  given  in  Fig.  14. 

4  Potential  Use  by  the  Federal  Government 

The  most  immediate  use  within  government-sponsored  programs  is  within  the  field  of  un¬ 
steady  flow  control.  The  instantaneous  sensitivity  can  be  used  to  modify  flow  features 
advantageously.  Immediate  applications  could  be  dynamic  flutter  analysis,  tactical  missile 
propulsion  and  control,  and  high-performance  aircraft  maneuverability. 
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Density  at  time  t  =  0. 1 42  Velocity  at  time  t  =  0. 142 


Figure  10:  Riemann  Problem,  Flow  Variables 


Density  Sens,  at  t  *  0.142  Velocity  Sens,  at  t  =  0.142 


Figure  11:  Riemann  Problem,  Sensitivity  Variables 
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Density  Sensitivity  Comparisons  Momentum  Sensitivity  Comparisons 


Energy  Sensitivity  Comparisons 


Figure  12:  Comparison  of  CSE  with  Central  Differences 


Figure  13:  Flow  and  Sensitivity  Solutions  for  Initial  Density  Sensitivity. 


Density  al  lime  l  *  0  142 


Figure  14:  Flow 


and  Sensitivity  Solutions  for  Bezier  Parameter  Sensitivity. 
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